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P-h \ Abstract 

In the relativistic Quantum Molecular Dynamics (RQMD) approach 
baryons are described within the framework of covariant hamilton con- 
straint dynamics. The inclusion of a relativistic mean field results in a 
quasiparticle picture for the baryons. This requires to distinguish be- 
£NJ , tween canonical and kinetic variables of the particles. As resonances we 

include the A(1232) and the iV*(i440) resonance. The resonance masses 
f ^ ■ are distributed according Breit-Wigner functions. However, the scalar 

t^- self energy leads to a shift in the masses and introduces an additional 

medium dependence. Consequences of this description on resonance and 
pion dynamics are discussed. 



1 Introduction 



In relativistic heavy ion collisions, e.g. in the SIS energy range excited nucleon 
states play a decisive role in the reaction dynamics. Most important is the 
A (1232) resonance which can reach abundances comparable to nuclear matter 
saturation density |l|, however also the iV*(1440) plays a non-negligible role. 
The decay and rescattering of these resonances are the predominant sources 
of meson production. In order to use such mesons as a source of information 
of the hot and compressed phase in heavy ion reactions the understanding of 
their origin, i.e. of the properties of excited nucleon states appear to be indis- 
pensable. Microscopic studies of the excitation and deexcitation of nucleons 
can be performed within the framework of transport models like Quantum 
Molecular Dynamics (QMD) and the BUU-type models. 

In the covariant extension of QMD, i.e., the Relativistic QMD (RQMD) 
[0, [| up to now only static Skyrme forces were used. These generalized Skyrme 
forces are treated as scalar potentials in the framework of Constrained Hamil- 



ton Dynamics. But the full Lorentz-structure of the nucleon-nucleon interac- 
tion contains large scalar and vector components. Effects of these relativistic 
forces were already studied in the framework of covariant generalizations of 
models of the BUU-type ]5], || [7|, mostly in the framework of the Walecka 



model || and its non-linear extensions, but not yet in RQMD. In Ref. [0| the 
formulation of a sclar- vector RQMD has been given. Here we study in partic- 
ular the influence of the relativistic mean field on the resonance dynamics. 

2 RQMD with scalar-vector mean fields 

The relativistic self-energy S = E s — 7 M S M contains scalar and vector com- 
ponents. Thus, one has to distinguish between canonical momenta pi and 
bare masses Mi on the one hand and kinetic momenta p* and effective (Dirac) 
masses m* on the other hand. The latter correspond to the quasiparticles 
dressed by the surrounding medium which obey the in-medium Dirac equation 

W-<K = o (i) 

where u* is an in-medium spinor or a Rarita spinor in the case of a A. 



In the formalism of Constrained Hamiltonian Dynamics ||, [Uj, the 8N 
dimensional phase space of N interacting relativistic particles is reduced to 
6N dimensions by 2N constraints which fix the individual energies and by the 
mass-shell constraints parametrize the world lines by fixing the relative time 
coordinates. 

With respect to Eq. (|]) the N on-mass-shell conditions for the four- 
momenta are given in terms of the quasiparticles 

^ = (rf-sr)(^-s i ,)-(M l -s«) 2 =^r-< 2 = o , (2) 

which contain scalar and vector self-energies. Hence, the constraints, Eq.(0), 
take the full Lorentz structure of the relativistic mean field into account. In 
the relativistic Hartree approximation the scalar and vector self-energies are 
proportional to the scalar density p s and the baryon current B 11 , respectively 

E s = T aPa , TP = T V B» . (3) 

The proportionality factors r Sj „ in Eq. (|j) may either be constants as in the 
Walecka model || or, if one includes higher order medium effects, density 
dependent functions as discussed, e.g, in Ref. fLlf . Here we apply the non- 
linear Walecka model (NL2) which can be expressed in the form of Eq. (|3J) 
with 

gl 

Ts{ps) = ml + B$(p°) + C&(p a ) (4) 



where $ is the scalar a-meson field [II]. The scalar density and the baryonic 
current 

Ps(qi,Pi) =Eft ' B MiiPi) =Y,PiJ u m ( 5 ) 

are determined in terms of the scalar two-body densities 

Pij = (4 7r L)3/2 eXp( ^ /4L) ' (6) 

with u^ = p^/m* being the 4-velocity of particle j and qxij the invariant 
center-of-mass distance of the particels i and j as defined in Ref . [[| . 

The total Hamiltonian is given by the N on-shell constraints, Eq. (0), and 
by N-l time fixation constraints </>, choosen as in ||. These time constraints 
ensure that interacting particles have equal times in their center-of-mass frame. 
A final constraint fixes the global time evolution parameter. The Hamiltonian 

N 2N-1 

H = "£\ i K i + Yl ^ (7) 

i=l i=N+l 

generates the equations of motions for the canonical conjugated coordinates 
and momenta 

dtf/dT=[H,qJt] , drf/dT = [H,tf] (8) 

where [A, B] means the Poisson bracket of phase space functions A and B. 
Notice that in the present formalism canonical and kinetic momenta are no 
independent quantities. To integrate the set of above equations @, one has 
to determine the unknown Lagrange multipliers \{r). This can be done using 
the fact that the complete set of 2N constraints must be fulfilled during the 
whole time evolution, i.e. dKi/dr = dfy/dr = 0. If the Dirac's first class 
condition is fulfilled, i.e. [Ki, Kj] = 0, the Hamiltonian is reduced to 

JV 

H = J2X i K i , Ai = A^ (9) 

where A^ = [Kj , Xi] is a submatrix of the complete constraint matrix Cji [BL ffl . 
In this approach the dynamics are to high extent determined by the La- 
grange multipliers A$(t). Thus in Fig.l we show the time evolution of A for a 
representative nucleon in a central Ca+Ca collision at 1.85 A.GeV for both, a 
pure mean field calculation and including binary collisions. It is clearly seen 
that A is a smooth function in time if its evolution is governed only by the mean 



field dynamics. However, in the high density phase (10/m/c < t < 20 fm/c) 
strong fluctuations occure which are due to binary collisions. These fluctua- 
tions represent rapid changes in the momenta and - in the case of a resonance 
creation or decay - the masses of the particles. 
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Fig. 1. Time evolution of a lagrange multiplier in a central Ca+Ca reaction at 
1.85 A.GeV. A full calculation including binary collisions (solid) is compared 
to a pure mean field calculation (dashed). 

The inclusion of a relativistic mean field which introduces an additional mo- 
mentum dependence into the model is essential in order to retain the corrcet 
dynamics at incident energies above about 1 A.GeV. In Fig. 2 we compare 
the present approach to the standard RQMD || |j] which uses static Skyrme 
forces. From Fig. 2 which shows the transverse flow per particle for the reaction 
Ar+KCL at 1.8 A.GeV under minimum bias condition it is clearly seen that 
the Skyrme forces result in far too less repulsive mean fields whereas NL2 is 
in good agreement with the data. 



3 Baryonic resonances 

The mean field propagation described in the previous section is the same for 
nucleons and baryonic resonances. In the latter case the bare nucleon mass in 
Eq. (Q) has to be replaced by the corresponding resonance mass Mr, i.e. m* R = 
Mr — S s r. Thus we apply the same mean field to nucleons and resonances, 
i.e., the coupling strenght of the respective mesons is assumed to be identical. 



Since we restrict to non-strange mesons this appears to be reasonable In the 
present approach the coupling functions T s >v have to be interpreted as effective 
quantities which parametrize the mean field rather than elementary vertex 
functions which further justifies this assumption. 
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Fig. 2. Transverse flow per particle for the reaction Ar+KCl at 1.8 A.GeV. 
Static Skyrme forces (soft/hard) are compared to the non-linear Walecka 
model NL2 and to data from Ref. [[12]. 

For the inelastic nucleon-nucleon channels we include the A(1232) as well as 
the A/"* (1440) resonance with the cross sections of Ref. [[H|. The lifetimes of 
the resonances are determined through their energy and momentum dependent 
decay widths 

ai|p|° 



(l + a 2 |p| 2 )(a 3 + |P| ) 
which originates from the p-wave representation of the resonances. In Eq. (|T0"D 
p is the momentum of the created pion (in GeV/c) in the resonance rest frame. 
According to Ref. |3| the values ai=22.83 (28.8), a 2 =39.7 and a 3 =0.04 (0.09) 
are used for the A (JV*) and the bare decay widths are taken as T^= 120 MeV 
and T^*= 200 MeV. 

Since only the quasiparticles, i.e. kinetic momenta and effective mass lie on 
the mass-shell all collisions are performed in the kinetic center-of-mass frame 

with a/s* = JpcM + m i 2 + v Pcm + m 2 2 - m the case of an inelastic collision, 
e.g. AW i — ► NR the final momentum can be evaluated 
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with A = [s* — (m* N + m R ) 2 ][s* — (m* N — m* R ) 2 }. The probability distribution 
for the effective resonance mass is given by a Breit-Wigner distribution 
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(12) 



with m° R = M R — S s r. The maximal possible bare mass of the created reso- 
nance is then restricted by 



M™ 



y/s* — m* N + S 
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sR 



(13) 
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Fig. 3. Shift of the maximal resonance mass by the scalar self-energy in nuclear 
matter (solid, left scale). In addition the effective nucleon mass m* is shown 
(dashed, right scale). The non- linear Walecka model NL2 is applied. 



Fig. 3 shows the medium dependence of M R ax created in a A^A^ i — > NR colli- 
sion. Here we have choosen the momenta of the incident nuclei (p* = — p£) to 
be 0.6 GeV/c in the nuclear matter rest frame. It is seen that M R ax is consid- 
erably enhanced by the presence of the medium, i.e. by the attractive scalar 
self-energies which is correlated to the reduction of the effective mass m* also 
shown in Fig. 3. Thus the probability for the excitation of higher resonances 
is generally enhanced in the presence of relativistic mean fields. This effect is 
even more pronounced when models are applied which result in larger fields 
than NL2 as, e.g., the original Walecka model M. 



The reabsorption processes (7rN — > A, N*) are treated as described in Ref. 



15| again adopting the cross sections from Ref. [JT3 

Finally in Fig.4 we show the pion p t -spectrum obtained in a 1.85 A.GeV 
Ni+Ni collision. The contributions from pions originating from iV*-resonances 
are shown seperately and it is seen that these pions contribute over all by about 
10% to the total yield. However, N* pions are in particular relevant for the 
high energetic part of the spectrum. Here the contributions from N* are more 
than 20%. Since high energetic pions are supposed to be most likely produced 
in the early phase of the reaction [|16J they really probe the high density phase 
and thus can give signals on the medium dependence of the resonances. 
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Fig. 4. Transverse 7r° spectrum in a Ni+Ni collision at 1.85 A.GeV under 
minimal bias condition. The contribution from pions originating from N* are 
shown seperately. 



To summarize we have applied relativistic mean fields with scalar and vector 
components in the formalism of Hamilton Contraint Dynamics. This results in 
a quasiparticle picture for baryons and resonances. The attractive scalar fields 
lead thereby to a shift of the resonance mass distributions towards higher val- 
ues and thus reduces the thresholds for the excitation of high lying resonances. 
Concerning heavy ion reactions in the SIS domain we found that the N* gives 
important contributions to the pion spectra. In particular high energetic pi- 



ons will be well adopted to study the influence of medium effects on baryonic 
resonances. 
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